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The validity of the adiabatic approach to spin-Peierls 
transitions is assessed. An alternative approach is devel- 
oped which maps the initial magneto-elastic problem to an 
effective magnetic problem only. Thus the equivalence of 
magneto-elastic solitons and magnetic spinons is shown. No 
soft phonon is required for the transition. Temperature de- 
pendent couplings are predicted in accordance with the ana- 
lysis of experimental data. 
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Around 1980 there was a great interest in the phe- 
nomenon of spin-Peierls (SP) transitions |1J where the 
coupling of lattice degrees of freedom to quasi one- 
dimensional [d = 1) magnetic degrees of freedom leads to 
a phase transition into a dimerized phase. This interest 
has been vividly renewed recently due to the discovery of 
the first inorganic SP substance CuGeC>3 ||||. 

The instability of the coupled spin-phonon system to- 
wards dimerization results from the susceptibility of the 
magnetic d = 1 subsystem towards 2fcp perturbations, 
i.e. dimerization. The gain in magnetic energy due to 
dimerization overcompensates the loss in elastic energy 
of the lattice distortion. The present work is motivated 
by the wealth of information available for the SP transi- 
tion in CuGeC>3. 

The theoretical picture of SP transitions has been de- 
velopped in a number of articles (e.g. |fl-^fj|). It relies 
so far mostly on an adiabatic treatment of the phonons. 
Cross and Fisher discussed this point most comprehen- 
sively (6). By random phase approximation (RPA) they 
investigate the stability of the uniform phase. The RPA 
is the consistent extension of the mean-field treatment on 
the one-particle level to the two-particle level of suscep- 
tibilities. Cross and Fisher's point is the importance of a 
11 pre- existing soft phonon" . This means that the phonons 
responsible for the distortion have a low energy already 
before the interaction with the correlated system is taken 
into account. Moreover, they require that the phonon 
dispersion perpendicular to the chain direction is very 
large. Thus the lattice prefers that whole planes perpen- 
dicular to the correlated chains move coherently so that 
the moving objects are heavy. 

The reason for the above requirements is that the 
mean-field approach is appropriate if the fluctuations are 
small compared to the expectation value. This is the 
case if the distortions are made up by a large number of 
phonons which in turn means that the phonon energies 
must be small. The phonons must be slow and heavy. 
Then it is plausible to consider the phonon subsystem 



as the slow subsystem which is renormalized by the fast 
magnetic degrees of freedom. In this picture, the SP 
transition is signalled by the vanishing of a renormalized 
phonon frequency H. 

Whereas the experimental data for organic SP sub- 
stances such as (TTF)CuS 4 C4(CF 3 )4 supports the RPA 
approach J(| , the experimental evidence for the inorganic 
CuGe03 points into the opposite direction |0|. Braden 
et al. found that two out of four optical phonons allowed 
by symmetry contribute appreciably to the distortion 
in a ratio of 3:2 JTTf . The more important phonon is 
at 6.8Thz (330K); its dispersion at the zone boundary 
is essentially flat decreasing towards the zone centre to 
w 3THz. The other phonon is at 3.2THz (150K) and 
practically dispersionsless. The magnetic exchange cou- 
pling J lies in the range 115K to 160K @-|||. It is ev- 
ident that the assumption of a pre-existing soft phonon 
is inadequate. Things get even worse if one takes into 
consideration the results for a XY chain with d = 1, 
Einstein phonons JHJ which show that the phonon en- 
ergy u> must be small compared to the resulting gap 
A for the adiabatic approach to be reasonable. The 
condition lo < A, however, is definitely not fulfilled in 
CuGe0 3 with A = 23K [p|. In view of these facts it 
is not astounding that so far no phonon softening at 
the SP transition was found experimentally. For com- 
parison, we recall the numbers for the best-known sub- 
stance (TTF)CuS 4 C4(CF 3 )4 with soft phonon: J = 68K, 
T SP = 12K (=*> A w 1.77T S p ~ 20K), to w 10K where 
lu < A is fulfilled §. 

In the present work we develop a route to Peierls tran- 
sitions not based on the assumption of phononic adia- 
baticity. Results of previous works |ll]-|2l| on nonadi- 
abaticity will be extended. We will view the phonon 
subsystem as fast and the spin subsystem as slow. The 
unperturbed spin system is always gapless, i.e. the mag- 
netic subsystem always has low-lying excitations well be- 
low any (optical) phonon energy. These low-lying excita- 
tions are influenced most by the interaction of phonons 
and spins. To them the phonons are fast. So we treat 
the phonons as quickly adapting and derive an effective 
dressed spin model. Pytte did the same for an Ising 
model which allowed the rigorous elimination of phonons 
|l9f stressing already the importance of avoiding a mean- 
field approximation for the displacements. 

Technically we use the recently developed flow equa- 
tion approach to treat the spin-phonon system p2]| . The 
idea is to rotate away the direct interaction with phonons 
similar to what is done in Frohlich's approach p3[ . This 
approach has been improved considerably by Lenz and 
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Wegner for electron-phonon interactions |24(| . In the im- 
proved version the generated effective couplings are by 
far less singular than in Frohlich's approach. 

The flow equation approach "diagonalizes" a Hamil- 
tonian in a continuous unitary transformation 
parametrized by / € [0,oo]. This means H(0) is the bare 
Hamiltonian as given and H(oo) is the resulting (more) 
diagonal Hamiltonian. The unitary transformation is de- 
fined by its antihermitean infinitesimal generator r](l) via 



dH 

A good choice for r) is 



[r)(l),H(l)] . 



r,= [Hv(l),H(l)] 



(1) 



(2) 



where Hd is the suitably chosen diagonal Hamiltonian 
p2| . The main feature of r/ as defined in (j^) is that it 
respects the idea of renormalization in that it eliminates 
first matrix elements connecting very different energies 
II- 

The SP system we consider is given by the Hamiltonian 



H = H S + H B +H SB H S =J2 d(q)S q ~S^ (3a) 
ff B = 5>W 6 ? H SB = J2 A M + b_^ (3b) 
A ? =J2Mk)S~ k S_~ k _ ? . (3c) 



in obvious notation in momentum space. Note that 
according to Ag should be normal-ordered Ag — > 
Aq- — {Ag) . The particular choice (before eq. |l2| and 
cq. [l4]) for Aq- will circumvent this problem. We as- 
sume inversion symmetry so that uj(q) = Lu(—q) and 
d(q) = d(—q). Hermiticity requires A~^ g — A- or equiv- 

alently g*{—q,—k) — g(q,k). The linear boson terms 
become /-dependent for the unitary transformation 



h sb (1) = J2(t ? (i)bi + T+(i)b cJ 



(4) 



with the starting condition T^-(O) = Ag. The objective of 
the unitary transformation is to disentangle phonons and 
spins. Thus we choose Hd = H§ + H B . To leading order 
in g we do not need to consider a possible /-dependence of 
-Hd since the /-dependent terms enter only in order g 2 as 
we will see. These induced terms of order g 2 and higher 
lead to a new contribution AH to the Hamiltonian. In- 
troducing the Liouville operator C for the commutation 
with Hq- CA := [iJs>^4] we choose for the generator rj 

V (l) = [H B ,H SB ] (5a) 
= ((£ + <"®)T f W + OC - "fflT+Wbf ) (5b) 



which is motivated by (@). The flow equation (|j) leads 
to 

flJ-f 

— = [ V ,H SB } + [ V ,AH} 

-£((£ + (Z)6t + (£ - o;(g)) 2 T+(/)& ? ) (6) 



^- = -{L + ^q)fT 3 



In linear order in g we have to fulfil the flow equation 

which is formally solved by 

T ff (l) = e X p(-(£ + u(q)) 2 l)A ff . (8) 

Based on (||) the additional Hamilton part AH can be 
calculated 



dl 



= [ V ,H SB ] + 0(g 3 ) 
^ |_ g <? ' 9 g ' 9 « 99 

q. k 

V (.D+T- + T+.D- 

l 9 9 9 9 



E(^[^.^] + h.c. 



(9a) 

G(ff 3 ) (9b) 
(9c) 
(9d) 



h th (PV + [^»^])J +°(5 3 ) (9e) 



q . k 



where we used the shorthand Dg = (C + ui(q))Tff. 

To obtain from (Q) an effective spin Hamiltonian we 
use a mean-field approach and replace the quadratic bo- 
son terms by its expectation values. This is absolutely 
systematic in the sense of an expansion in g. Taking the 
expectation values neglects fluctuation effects of the or- 
der g 2 due to the interaction. But since the two-boson 
terms appear only as g 2 -terms the total error due to the 
mean- field treatment is of the order g 4 . Applying the 
same mean-field approach to the unspecified g 3 terms an- 
nihilates them because they contain necessarily an odd 
number of boson operators. Thus the effective spin model 
is exact up to 0(g 4 ). 

Replacing btb- by ^(ex-p(w(q)/T) — 1) _1 and omit- 
ting the terms in ( |9c| ) we obtain finally 

dAH ^AV+coth(^iV^ with (10a) 



dl 



Q 

1 



2T 



(10c) 



^=-^(P 9 T+ + D+T 9 (10b) 

Yo = l([T ff ,Di] + [D q -,T+] 

From (|h]) the Hamiltonian corrections AH x /y are found 
by integration over / and summation over q: AHx = 
Jo°° Y, 3 X3dl and AH Y = /~ ^-.coth (u(q)/ (2T)) Ygdl. 
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In order to get an impression of what @ means we 
assume J<u(q) and calculate the leading contributions 
in C to X$ (even in C) and Y$ (odd in £). After some 
algebra and integration we find 



AH X = - V-J-A+A. 



„ ; — coth , 



uj(q) 



(lla) 
(lib) 
(11c) 



where we simplified the formulae (Jl la|) and (11c) in real 



space one step further approximating the phonons by 
Einstein phonons. The term AHx corresponds to the 
results obtained previously by other methods |l^ |2l|] . To 
the author's knowledge, the T-dependent term AHy has 
not yet been described. The result of Pytte []l9| is found 
back by observing that in the Ising model CAi vanishes 
since all terms involving only S z commute, i.e. AHy be- 
comes zero. The result in pl| | is retrieved on observing 
that g 2 /uj is proportional to J 2 /(muj 2 ) in |0j since the 
displacments equal Ui = (b{ + b])/^/2muj. Neglecting the 
phononic kinetic energy while keeping their potential one 
constant corresponds to the limit m — ► with raw 2 con- 
stant, i.e. u> — > oo. So the g 2 /ui term is constant and kept 
while terms g 2 0{uj~ 2 ) like AHy are neglected in ^J. To 
further enhance the plausibility of the result ([ll]) we note 
that it equals the result one gets by Frohlich's method 
p3[ in the two leading orders g 2 /oj and g 2 J/uj 2 . The 
difference of the flow equation approach and Frohlich's 
approach appears only in the 1/w 3 terms coming from 
©. Insofar Frohlich's approach can also be used to 
derive ([ll]). The flow equation approach, however, is a 
better starting point for future higher-order calculations 
in g/to and J/ui which take ^-dependent couplings into 
account. This is the reason why this method is chosen 
here. 

Specifically, we consider first strictly one dimensional 
phonons A4 = g(SiSi+i — SiSi-±). This choice guarantees 
(Ai) = in the symmetry unbroken phase so that is 
normal-ordered. With this A- we have 



AHy 



9 2 ^ 



{S{Si- 



- Si S; 



(12) 



For AHy we have to know d(q) in ([3a]). We assume 
nearest and next-nearest neighbour interaction J and aJ, 
respectively: d(q) = J(cos(gi) + acos(2gi)). One obtains 



AHy = — --r coth ( — ) x 



4 uj 2 \2T) 
^(-(3 - 3a)SiS l+1 + (3 - 5a)SiS i+ 2 + 2aS l S i+3 ) 

i 

where products with four different spins are omitted. 



(13) 



Even if no frustration is present in the original model 
(a = 0) the dressing of the spins with phonons induces 
frustration a c s > 0. The couplings are temperature de- 
pendent since they are mediated by the phonons. Using 
the term "spinon" for a purely magnetic elementary S — 
1/2 excitation and the term "soliton" for the joint mag- 
netic and elastic S = 1/2 excitation we state that the soli- 
tons of the Hamiltonian (^) are unitarily equivalent to the 
spinons of the Hamiltonian H = Hs+Hb + AHx + AHy. 
This shows that solitons and spinons are in essence the 
same entity and puts Affleck's supposition in this respect 
p6| on a quantitative basis. 

The low-lying excitations of the frustrated Heisenberg 
chain are spinons which are gapless for a c ff < a c = 0.241 
p7f and gapful above a c |28 where the system undergoes 
also a spontaneous symmetry breaking of the transla- 
tional symmetry towards a dimerized phase. The contin- 
uum starts right at the gap energy Q . These facts imply 
that a single chain shows a SP transition only above a 
certain value of the interaction in contrast to the results 
of the adiabatic treatment (see also |]l7| ). Furthermore, 
no "double gap" feature |3C| occurs in a single chain. 

An elastic interchain coupling in a chain ensemble is 
included if the local phonons influence also neighbouring 
chains 

^-i,3 = 9\&i,j ~ Si,jSi—i,j + 
<j,3'> 

where j is the chain index and j and j' are adjacent 
chains. The factor |/| < 1 indicates the influence of 
a certain distortion on one cha in onto adjacent chains. 
Due to the commutators in (11c) a finite / changes AHy 

5 2 (1 + zj 2 ) 



<r 



in (y_3|) only by renormalizing g 
where z is the number of neighbouring chains each chain 
has. The same renormalization takes place in AHx- 
Additionally, terms linking different chains occur like 



'{9 f /^{SijSi+i j — SiSi-±)(Si t j> Si-ij' 



■1,J' 



These terms drive really the finite temperature SP tran- 
sition since they enable at low enough temperature a co- 
herent dimerization throughout the whole lattice. We 
will call these terms coherence terms. Their influence 
on the low-lying excitations is to confine pairs of spinons 
(solitons) to triplets or to singlets. For the realistic case 
of weak coupling g 2 f '/u> <C J a mean-field treatment is 
justified. This amounts up to the treatment of dimer- 
ized, frustrated chains with self-consistently determined 
dimerization. Hence the confinement is the same as in 
dimerized chains (see e.g. 26 3J]). This explains why 



the adiabatic approaches based on dimerized, frustrated 
chains capture correctly the physics of the dimerized SP 
phase at low T. 

The main difference of the fast-phonon scenario to the 
adiabatic one is the absence of a soft phonon at the tran- 
sition. The transition is characterized by growing do- 
mains of coherent dimerization the size of which diverges 
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at Tsp. No renormalized phonon frequency needs to van- 
ish. Interestingly, the RPA shows similar results in the 
nonadiabatic parameter regime Q . 
How does this fit to the approach used so far? In the 
usual RPA treatment j(| the phononic self-energy con- 
tributes not only a real part but also an equally strong 
imaginary part which stands for strong damping p3l . 
Thus the real, untransformed phonons are not appropri- 
ate quasi-particles. If we transform the propagator of the 
real phonons in the same unitary way as the Hamiltonian 
we see that via b\ — > b\ — ^Ai + . . . not only the trans- 
formed phonons matter but also S = excitations of the 
effective spin model. This means that the observation 
of the real phonons reveals not only a sharp peak at the 
high frequency to but also a continuum at low frequencies 
of the order of J. The low energy continuum changes on 
passing through the spin-Peierls transition. In the vicin- 
ity of Tsp one expects some critical fluctuations close to 
zero energy. Below Tsp a gap should appear which equals 
twice the triplet gap or less if bound states are present 
|ii"| , j34| . But the phonon peak is not lowered towards zero 
energy at the transition. This is in accordance with the 
results known so far for CuGeC>3. 

We attempt to estimate orders of magnitude for the 
couplings in CuGeC>3. Let us assume that Tsp is of the 
order of g 2 /uj then g 2 /u 2 w T SP /u w 15K/150K = 0.1 is 
roughly one tenth which justifies the expansion in g 2 /uj 2 . 

The extension of ([ll]) to several phonons is straight- 
forward. Using the phonon energies and their relative 
distortion as experimental input and the coupling g as 
fit parameter it is possible to reproduce the experimen- 
tal x{T) data nearly as well as in jl6| from where also the 
x(T) data was taken. This shows that the assumption 
of T-dependent couplings is not ruled out by the x{T) 
data. With the same parameter as for the x(T) fit we find 
J(50K) = 162X and J(300K) = UOK which agrees very 
well to J(50K) = 158X and J(300K) = 136X deduced by 
Fabricius and Low f35f| from experimental S(q,u) data. 
This excellent agreement confirms the validity of the ap- 
proach used and in particular the prediction of T depen- 
dent couplings. 

In summary, we discussed the validity of the phonon 
adiabatic approach for SP transitions and in particu- 
lar CuGeOs. The phonon adiabatic approach is inad- 
equate for the latter system. We developed a promising 
alternative approach relying on the flow equation tech- 
nique. Magneto-elastic solitons are mapped to magnetic 
spinons in an effective, purely magnetic Hamiltonian. 
The phonon dynamics induces a T-dependent frustration. 
No soft phonon signals the SP transition which is driven 
by the coherence terms in an effective magnetic model. 

The author is indebted to M. Braden for making his 
results available to him prior to publication. He acknowl- 
edges helpful discussions with M. Braden, W. Brenig, B. 
Biichner, E. Muller-Hartmann, and H. J. Schulz. This 
work was supported by the DFG through SFB 341. 



[1] J. W. Bray, L. V. Interrante, I. C. Jacobs, and J. C. Bon- 
ner, in Extended Linear Chain Compounds, edited by 
J. S. Miller (Plenum Press, New York, 1983), Vol. 3, p. 
353. 

[2] M. Hase, I. Terasaki, and K. Uchinokura, Phys. Rev. 

Lett. 70, 3651 (1993) 
[3] J. P. Boucher and L. P. Regnault, J. Phys. I France 6, 

1939 (1996) 
[4] E. Pytte, Phys. Rev. B 10, 4637 (1974) 
[5] J. W. Bray et ai, Phys. Rev. Lett. 35, 744 (1975) 
[6] M. C. Cross and D. S. Fisher, Phys. Rev. B 19, 402 

(1979) 

[7] M. C. Cross, Phys. Rev. B 20, 4606 (1979) 

[8] T. Nakano and H. Fukuyama, J. Phys. Soc. Jpn. 49, 1679 

(1980); ibid. 50, 2489 (1981) 
[9] A. I. Buzdin, M. L. Kulic, and V. V. Tugushev, Solid 
State Commun. 48, 483 (1983) 
[10] M. Fujita and K. Machida, J. Phys. Soc. Jpn. 53, 4395 

(1984); ibid. J. Phys. C21, 5813 (1988) 
[11] M. Braden et al, Phys. Rev. B 54, 1105 (1996); ibid. 

Phys. Rev. Lett, in press 
[12] M. Nishi, O. Fujita, and J. Akimitsu, Phys. Rev. B 50, 
6508 (1994) 

[13] J. Riera and A. Dobry, Phys. Rev. B 51, 16098 (1995) 
[14] G. Castilla, S. Chakravarty, and V. J. Emery, Phys. Rev. 

Lett. 75, 1823 (1995) 
[15] G. S. Uhrig, Phys. Rev. Lett. 79, 163 (1997) 
[16] K. Fabricius et al, Phys. Rev. B 57, 1102 (1998) 
[17] L. G. Caron and S. Moukouri, Phys. Rev. Lett. 76, 4050 

(1996) 

[18] L. P. Regnault et al, Phys. Rev. B 53, 5579 (1996) 
[19] E. Pytte, Phys. Rev. B 10, 2039 (1974) 
[20] U. Brandt and H. Leschke, Z. Phys. 271, 295 (1974) 
[21] F. H. L. Essler, A. M. Tsvelik, and G. Delfino, Phys. Rev. 

B 56, 11001 (1998) 
[22] F. J. Wegner, Ann. Physik 3, 77 (1994) 
[23] H. Frohlich, Phys. Roy. Soc. Lond. A215, 291 (1952) 
[24] P. Lenz and F. Wegner, Nucl. Phys. B 482 693 (1996) 
[25] S. K. Kehrein and A. Mielke, Ann. Phys. 252, 1 (1996) 
[26] I. Affleck, in Dynamical Properties of Unconventional 

Magnetic Systems (NATO ASI, Geilo, Norway, 1997). 
[27] R. Jullien and F. D. M. Haldane, Bull. Am. Phys. Soc. 

28, 344 (1983); K. Okamoto and K. Nomura, Phys. Lett. 

A169, 433 (1992) 
[28] R. Chitra et al, Phys. Rev. B 52, 6581 (1995); 

S. R. White and I. Affleck, Phys. Rev. B 54, 9862 (1996) 
[29] B. S. Shastry and B. Sutherland, Phys. Rev. Lett. 47, 

964 (1981) 

[30] M. Am et al, Phys. Rev. Lett. 78, 1560 (1997) 
[31] G. S. Uhrig and H. J. Schulz, Phys. Rev. B 54, R9624 
(1996) 

[32] H. J. Schulz, private communication 1998 

[33] C. Gros, Korrelationstage, Dresden 1997 

[34] G. Bouzerar, A. P. Kampf, and F. Schonfeld, 

cond/mat9701176 
[35] K. Fabricius and U. Low, |sond-mat 9709160 



4 



